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Y2 ' Abstract 
(N 

\ rea, where the off-axis neutrino beam from J-PARC at Tokai village for the T2K 

' project has significant intensity at a few GeV range. In particular, we examine 

the capability of determining the mass hierarchy pattern and the CP phase of the 
\ lepton-flavor-mixing matrix when a 100 kt water Cerenkov detector is placed at 

I ' various locations in Korea for the off-axis beam (OAB) of 2.5° and 3.0° at the 

Super-Kamiokande site. The best results are found for a combination of 3.0° OAB 
r~| ! at SK (L = 295km) and 0.5° OAB at L = 1000km, where the mass hierarchy 

pattern can be determined at 3-0" level for sin^ 20rct ~ 0.05 (0.06) when the hier- 
archy is normal (inverted), after 5 years of running (5 x 10^^ POT). We also find 
that the leptonic CP phase, Jmns , can be constrained uniquely, without invoking 
I anti-neutrino beams, as long as the mass hierarchy pattern is determined. Those 

results are obtained by assuming that the charged current quasi-elastic events can 
be separated from the other backgrounds with high efficiency, the neutrino en- 
ergy can be reconstructed with a hundred MeV uncertainty, and the earth matter 
density along the baseline can be determined with 3% accuracy. 
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We study physics potential of placing a far detector in the east coast of Ko- 
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1 Introduction 



All the experimental results of the neutrino flavor oscillation [1] are consistent with 3 
neutrinos except for the LSND [2] experiment. Among the 9 parameters of the three 
neutrino model, 6 parameters can be measured by neutrino oscillation experiments: 2 
mass-squared differences, 3 mixing angles and 1 CP violating phase. As of 2006 spring, we 
have already known 2 mass-squared differences and 2 mixing angles from the atmospheric 
neutrino experiments [3], the long-base-line (LBL) neutrino oscillation experiments [4,5], 
the solar neutrino experiments [6,7] , and the reactor neutrino experiments [8,9] . However, 
the sign of the larger mass squared-difference (m| — mf), one of the 3 mixing angles 
(^rct), and the CP phase (5mns) have not been measured yet. The tasks of the future 
neutrino oscillation experiments are not only to confirm the 3 neutrino model, but also 
to measure those unknown parameters of the model. 

Here we focus our attention on one of the future neutrino experiments, the Tokai- 
to-Kamioka experiment (T2K) [10] which will start in 2008. The center of the T2K 
neutrino beam from J- PARC [11] at Tokai village will go through underground beneath 
Super-Kamiokandc (SK), and reach the sea level east of Korean shore. At the baseline 
length L = 295 km away from Tokai village the upper side of the beam at 2° to 3° 
off-axis angle is observed at SK, and the lower side of the same beam at 0.5° to 3.0° 
off-axis angle can be observed in Korea [12]. An additional far detector in Korea can 
probe the neutrino oscillation at a baseline length (L) of 1000 to 1200 km away from 
the Tokai village [12-14]. The most welcome feature of this two detector system has 
been identified [14] as the relative hardness of the neutrino beam energy spectrum at 
a smaller off-axis angle, which can be observed in the east coast of Korea at L ~ 1000 
km. Accordingly, it is possible to arrange such that a far detector in Korea probe the 
oscillation at around the same oscillation phase, 5m'^L/2E, at Kamioka. The difference 
between the resulting two measurements should then come from the difference in the 
matter effects, which can be a factor of three larger in Korea. In Ref . [14] , we showed 
that the two detector system can resolve the mass hierarchy ambiguity by making use 
of the strong matter effects on the v^j, — > Vf, transition probability [15-18] if a 100 kt 
level water Ccrcnkov detector is placed in the east coast of Korea and if the third mixing 
angle is not too small. It has further been shown in [14] that the leptonic CP phase, 
^MNS ) can be uniquely constrained once the hierarchy pattern is determined. 

In this paper we present details of our findings in Ref. [14]. The paper is organized 
as follows. In section 2, we review the formalism of the neutrino oscillation including 
the matter effect, as well as the constraints in the model parameters from the present 
experimental results. In section 3, we briefly introduce the T2K experiment and discuss 
merits of the Tokai-to-Kamioka-and-Korea (T2KK) proposal, where an additional large 
detector is placed in Korea along the T2K neutrino beam baseline. In section 4, we 
explain details of our analysis method where we select the charged current quasi-elastic 
events, and a function is proposed that takes into account statistical and systematic 
uncertainties in the experiment, as well as a part of possible background contaminations. 
In section 5, we show the results of our numerical calculation on the determination of 
the neutrino mass hierarchy. In section 6, we present our studies on the determination 
of the CP phase. In section 7, we summarize our findings and give some discussions on 
the possibility of expanding the T2KK two detector system, and the necessity of further 
studies on backgrounds which we could not include in the present analysis. 
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2 Oscillation formulae under the experimental constraints 



The neutrino flavor eigenstates, |z/q,) (a = e, yU,r ), are related to the mass eigenstates 
[i = 1,2,3) through the lepton-flavor-mixing matrix, or the Maki-Nakagawa-Sakata 
(MNS) matrix [19] 

Wa) = Uai Wi) . (1) 

The probability that an initial flavor eigenstate z/q, with energy E is observed as a flavor 
eigenstate z/^ after traveling a distance L in the vacuum is expressed as 

= \Up,U:, + Uf,,U:,e-^^^^ + Up,U:,e-'^-^\' , (2) 

where the phase Aj^ is 

A.- = -^^L ^ 2.534 ^^^^^^^ L[km] . (3) 

Eq. ([2]) tells us that neutrino oscillation experiments measure the 2 mass-squared dif- 
ferences and the lepton-number conserving contributions of the MNS matrix elements, 
which can be parametrized by 3 mixing angles and 1 CP violating phase. 

The present neutrino oscillation experiments are each sensitive to only one of the 
2 mass-squared differences, and the amplitude of each oscillation probability can be 
expressed in terms of the MNS matrix elements. The atmospheric neutrino oscillation 
experiments [3] and the long-base line (LBL) neutrino oscillation experiments, K2K [4] 
and MINOS [5], which measure the z/^ survival probability, are sensitive to the magnitude 
of the larger mass-squared difference. The constraints on the mass-squared difference 
and the amplitude are [3-5] 

1.5 X lO-^eV^ <\ml-ml\< 3.4 x lO^W , (4a) 
sin^ 2^ATM > 0.92 , (4b) 

each at the 90% confidence level. The solar neutrino oscillation experiments [6,7] and the 
KamLAND [8] experiment, which measure the i^e and Pe survival probability, respectively, 
are sensitive to the smaller mass-squared difference. The present constraints can be 
expressed as [7] 

ml-ml = (8.0 ± 0.3) x IQ-^eV^ , (5a) 
sin^ ^soL = 0.30 ± 0.03 . (5b) 

In the solar neutrino experiments, the sign of — mf is determined by the matter effect 
in the sun [20,21]. The reactor Ue experiments at L ~ 1 km are sensitive to the oscillation 
with the larger mass-squared difference. No reduction of the survival probability has 
been observed, and the CHOOZ experiment gives the upper limit on the amplitude [9]; 



sin^ 26'rct < 0.20 for |m| - ml 
sin^ 26'rct < 0.16 for \m\ — m\ 
sin^ 26'rct < 0.14 for jmo - m\ 



2.0 X 10-W^ , 
2.5 X lO^W^ , 
3.0 X 10- W 



(6a) 
(6b) 
(6c) 
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at the 90% confidence level. 

These observed amplitudes can be identified with the MNS matrix elements as follows 
[22]. In the atmospheric neutrino experiments and the CHOOZ reactor experiment, 
IA13I ~ 1 ^ A12 is satisfied. Therefore we can set e"*^^^ ~ 1 in eq. ([2]) and we find the 
following relations; 

l^^^tsp = sin^ 6* ATM , 

|f/e3r = sin2^RCT. (7) 

As for the solar neutrino experiments and the KamLAND experiment, where the A12 ~ 1 
region is probed, the terms with A13 oscillate rapidly within the experimental resolution 
of L/E. After averaging out the A13 contribution, and by neglecting the term of order 
|f/e3p, which is constrained to be smaller than about 0.04 by the CHOOZ experiment [9], 
eq. fl6b|l . we obtain the relation; 

4|f/elf/e2r = sin2 2^SOL. (8) 

These simple identifications, eqs. ([7]) and ([HD, are found to give a reasonably good de- 
scription of the present data in dedicated studies [23] of the experimental constraints in 
the three neutrino model. In this paper we parametrize the MNS matrix in terms of 
the three positive numbers, sin^ 6'atm, sin^ 6'rct, and sin^ 26'sol with the identification of 
eqs. ([7]) and ([8]), respectively, and the CP phase 



mm 



- arg f/es . (9) 



This convention [22] allows us to express the MNS matrix in terms of the three observed 
amplitudes, directly. 

The neutrino oscillation through the earth is complicated by the fact that z/g and z/g 
have the extra W-hoson exchange interactions with electrons in the matter [20]. The 
effect is small at low energies, and an approximation of keeping only the first and second 
order corrections in the matter effect and the smaller mass-squared difference has been 
found useful for analyzing the LBL experiments at sub GeV to a few GeV range [24,25]. 
Although we evaluate the oscillation probabilities numerically in this report, the following 
analytic expressions are found to be useful for the T2KK setup where the earth matter 
effects can be treated pertubatively. For the survival probability, an accurate analytic 
expressions is found by retaining only those terms linear in A12 and the matter effect 
term [14,25]: 

P^.^u, = 1 - sin2 2^ATM (1 + A^^) sin^ + fi^) . (10) 

Here and i?^ are the corrections to the amplitude and the oscillation phase, respec- 
tively, which are linear in the smaller mass squared difference and the matter effect term. 
These can be expressed as, 

aL 1-2 sin^ 6'atm . 2 n 
A13E cos^ 6'atm 



^ = — -: ^ — sm 6't;_ 

/ 9 „ \ vr L /sin^2^RCT\ 

-0.005 (l-2sin^.ATM)^^(^^j 



:ila) 
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_ aL 1 - 2 sin' ^atm . 2 n 

4£/ COS^ C^ATM 

(cos' 6'soL + tan' 6'atm sin' 6'sol sin' 6'rct - tan 6'atm sin 26'sol sin 6'rct cos 5mns 



/sin'2eRCT\'^' e- 
0.037 - 0.008 ( ^ ^^""^ \ cos 5mns 



lA 



13 1 



TT 



(lib) 



Here the term a gives the contribution of the extra potential for v^. 



2V2GFEne ^ 7.56 x 10" W ( — ^ I ( t^) 

Vg/cm 



where Ue is the electron number density and p is the matter density. In the second 
lines of eqs. (11 lap , (lllbp . we employ the mean values of the atmospheric and the solar 



oscillation parameters in eqs. (jl]) and ([5]), and p = 3.0g/cm . The simple analytic 
expressions of eqs. (fTO!) and (fTTj) reproduces the survival probability accurately with 
1% error throughout the parameter range explored in this report, except where the 
probability vanishes. 

According to eqs. (Illal) . (l4bll . and ( l6bl) . the magnitude of around IA13I = n 
should be smaller than about 2 x 10"^ (8 x 10~^) at L = 295 km (1000 km), and 
hence the amplitude of the survival probability is not affected much by the matter 
effect. Therefore, we can measure sin' 26'atm rather uniquely from the z/^ disappearance 
probability independent of the neutrino mass hierarchy and the other unconstrained 
parameters. The phase-shift term B'^ affects the measurement of Irrig — m'|, whose 
preferred value grows (decreases) for the normal (inverted) hierarchy, and the magnitude 
of the shift can be about to 2% (3%) when cos(5mns = 1 (—1) for sin' 26'r,ct = 0.1. In 
other words, unless we determine the hierarchy, a few percent level of uncertainty should 
remain as a systematic error of jm' — m'|. 

For the z/^ Vf, transition probability, we need not only the linear terms of A12 and 
a but also their quadratic terms to obtain a good approximation. We find 

P,^^,^ = 4 sin' ^ATM sin' ^rct { (1 + A') sin' (^^) + 5^ sin A13 1 + C% (13) 
where A'^, B^, and are the correction terms: 

^e oQ siu 2^sOL . r f ttL 

A = - — - cos 26'rct —- — 7, ^ sm ^mns 1 + 



A13E 2 tan 6'atm sin 6'rct V 2Ai3£^ 

A12 / aL\ ( sin2gsoL . ■ 2 n \ 

H T U^12 + T a ^ cos ^MNS - 2 sm 6'sOL 

4 V 2hJ ytant^ATM smt^RCT / 

1 /aL\' 3 / aL \' ,^ ^ , 
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0.11 



TT 



Ai3 295km 



\ - ?ili3!£L] _ 0.014 f-^)' + 0.0087 f-^-^V 

2 ; V295km; VA13 295kmy 



-0.29 



smW^± 0.054^ 



sin^26'RCT/ 



+0.015 



oe oz) , ^12 / Sm^ysOL 

B' = -— cos 20RCT + — ; — ^ 



0.10 , 

cos ()mns - 0.11 



sin 26'rcT' 



TT 



295km y vr 



(14b) 



A12 / sin26'soL 



cos Suns - 2 sin^ (^1 + ^ 



aL \ sin 20sol 



2i?/ tan fc^ATM sm t^RCT A13 \2E J 



1 faiy 



I i A A _ _ 0.0091- ' I ^ 



-0.080 — 

V 295km/ 



Ai3 V295kmy 



(14c) 



+0.15 



0.10 



.sm^26'RCT 



cos(5mns - 0.11 



IA13I L \ 

^^±0.054; * 

TT 



295km J 



/ 0.10 

+0.0072 sin 5mns 

Vsm2 2^RCT/ V TT 



+ 1.7 



L \ |A 



\ 1^131 

295km ) TT 



^ - ^ ^ sin 26'soL sin 26'atm sin 6'rct cos 5mns 



_ ^12 -^2 cya ^^^2 /I ^12 al/ 

(-^ — — sm /lysoLCOs ^atm 7:;— 7n=; 

4 2 2i? 



^ ah \ ^ 

+ ( ^ j sin^ 6'rct sin^ 6'atm 



(14d) 
(14e) 



0.0011 



/ L \ ^ sin^ 26'rct 

+0.00036 — — . 

V295km/ 0.10 



(14f) 



Here ± in the first lines of eqs. ( 114b[) and ( 114dll follows the sign of A13. The first 
and second terms in eqs. f ll4ap and f ll4cp are linear, and the other terms including 
eq. (I14ep are quadratic in A12 and a corrections. The quadratic terms can dominate the 
probability when sin^ 2^rct , which gives the leading term, is very small. We find that 
the above analytic expressions, eqs. (fT3l) and (fT4|) . are useful throughout the parameter 
range of this report, down to sin^ 2^rct = 0, except near the oscillation minimum where 
the approximation can give negative probability. The largest deviation is found for 
sin^ 26'rct ~ 0.005, 5mns ~ 0° and the inverted hierarchy, where the analytic formula 
underestimates the probability by about 5% for L = 1000 km. 

The term proportional to sin^(-^) in eq. ([T3|) is the main oscillation term in our 
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approximation, and A'^ shifts the magnitude of the amphtude. From eq. fll4bp . we find 
that the amphtude of the z/^ z/g transition probabihty is sensitive to the mass hierarchy 
pattern, because the first term changes sign. Its magnitude increases (decreases) by about 
10% for the normal (inverted) hierarchy at L = 295 km for sin (5mns = 0. The difference 
between the two hierarchy cases grows with L when L/E is fixed at |Ai3| ~ vr, reaching 
about ±40% at L = 1000 km. The shift is also sensitive to sin^MNS? which can decrease 
(increase) the amplitude by as much as 30% for sin^MNS = 1(— 1) when sin^ 26'rct = 0.1, 
independent of the mass hierarchy. 

The term proportional to in eq. ( fT3l) can be re-organized as 

P,^^,, ^ 4 sin^ ^ATM sin^ O^ct (1 + A') sin^ + + , (15) 

just like the z/^ survival formula, eq. (fTOj) . when 1^4*^1 and |i?^|are small. Since the term 
B'^ shifts the oscillation phase in this limit, we call it as the phase-shift term. It shifts the 
oscillation peak energy higher (lower) for the normal (inverted) hierarchy, by about 5% 
at L = 295 km and by about 20% at L = 1000 km, for cos^mns = 0. The cos^mns 
dependence of the shift is also significant, which can be ±10% for cos^mns = ±1 
if sin^ 2^RCT = 0.1. Note, however, the factorized expression eq. f|T5|) seizes to be a 
good approximation when l^^l or \B'^\ gets larger than about 0.5, which happens when 
sin2 20RCT ~ 0.01; see eqs. ffTifel) and ffTidl) . 

The above observation inspires a two detector system, where both detectors can mea- 
sure z/^ — >■ Ue oscillations around the oscillation maximum (|Ai3| ~ vr) but at significantly 
different base-line lengths. If the magnitude of the product of sin^ ^atm sin^ ^rct is large 
enough that the z/^ z/g oscillation is observed at both detectors, then the difference of 
the observed oscillation probabilities at the two locations determines the mass hierarchy 
uniquely, and hence also the product sin^ 6'atm sin^ 6'rct rather independent of ^mns- 
Once the hierarchy is determined, the measurements of the magnitude and the phase of 
the z/^ z/g oscillation measure sin^MNS and cos^mns, respectively, and hence the CP 
violating phase 5mns can be determined uniquely. In the following section, we will find 
that the T2KK two detector system can indeed satisfy the above conditions, and our 
analysis based on the exact numerical evaluation of the oscillation probability supports 
the simple picture presented in this section based on the approximate formula eqs. flTUl) 

-m- 

3 Merits of detecting the T2K off-axis beam in Korea 

Many neutrino oscillation experiments [26-32] are planned to measure sin^ 2^rct • Tokai- 
to-Kamioka (T2K) neutrino-oscillation experiment which will start in 2008 is one of them. 
In T2K, high intensity z/^ beam is produced by the proton accelerator which is under 
construction at J-PARC in Tokai-village. These z/^'s are shot to Kamioka, 295 km west 
from Tokai. Super-Kamiokande (SK) will measure both the z/^ survival rate and 

the z/^ — > z/g transition rate. The measurement of the z/^ — > z/g transition probability is 
the main purpose of the T2K experiment, because it tells us the magnitude of sin^ ^rct, 
see eq. f|T3l) . the last unmeasured mixing angle of the 3x3 MNS matrix. In order to 
observe the signal of z/^ — > z/g clearly, the neutrino beam should satisfy the following 
conditions. 
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1. The neutrino energy near the oscillation maximum (|Ai3| ~ vr) at L = 295km is 
expected to be around 0.5 GeV to 0.7 GeV according to the present experimental 
bound, eq. (jl]). The flux should hence be large in this energy region. 

2. High energy neutrinos produce vr'^'s via neutral current, which become background 
to the Ve Charged Current Quasi-Elastic (CCQE) events. Therefore, the flux of 
beam should be small at high energies. 

T2K adopts the off-axis beam which satisfles the above requirements [10,33]. We show 
the flux of the T2K off-axis z/^ beam [33] in Fig. [1] (a), for lO^^ POT/yr at L = 295 
km for various off-axis angles between 0° and 3°. It is clearly seen that the flux peaks 
at 0.55 to 0.75 GeV at 2° to 3° off-axis angles. In Fig. [H^b) we show the cross section 
per nucleon of the z/g and CCQE events off the water target [4], and in Fig. [^c), we 
show the product of the Ve CCQE cross section and the flux at 295 km for various 
off-axis angles. Because the neutrino energy reconstruction is essential to determine the 
oscillation phase, we use only the CCQE events in our analysis. Fig. [T](b) shows that 
the CCQE cross sections grow quickly above the threshold, become ~ 3.5 x 10~^^ cm^ 
at Ey ~ 0.6 GeV, and stay approximately constant at ~ 4.5 x 10^ cm at El/ ^ 0.8 
GeV up to ~ 5 GeV where the flux diminishes. We also show the typical — z/g 
transition probability at L = 295 km and that at L = 1000 km in Fig.[T](d), calculated for 
ml-ml = 2.5 x IQ-^eV^, ml-m\ = 8.2 x 10-^eV^ sin^ 2^atm = 1-0, sin^ 2^sol = 0.83, 
sin^ 26'rct = 0.10, 5mns = 0°, and p = 2.8g/cm^ for L = 295 km, and p = 3.0g/cm^ for 
L = 1000 km. From Fig. [1] (a), (c), and (d), we conflrm that the 2.0° to 3.0° off-axis 
beam (OAB) has a strong flux peak where the oscillation maximum is expected at SK. 
We also note that the beam at smaller off-axis angles has signiflcant flux in the 1.4 ~ 3 
GeV region where the oscillation maximum of the z/^ v^. transition is expected at 
L = 1000 km. 

During the T2K experimental period, the center of the z/^ beam from J-PARC goes 
through Kamioka at 2.0° to 3.0° beneath SK, and the lower side of the same beam will 
appear in Korea at various off-axis angles. We show in Fig. [2] the off-axis angle of the 
beam at the sea level for the 2.5° off-axis beam (a) and 3.0° off-axis beam (b) at SK. The 
baseline lengths are shown by vertical contours, and the off-axis angles at the sea level are 
shown by elliptic contours between 0.5° and 3.0°. The SK is slightly off the corresponding 
contour because it is about 320 m above sea level. We flnd that, when the 2.5° (3.0°) 
off-axis beam reaches SK, the 1.0° (0.5°) off-axis beam appears in the east coast of 
Korea which is 1000 km away form J-PARC. If a huge neutrino detector is constructed 
in the east coast of Korea, along the T2K beam direction, the two detector system, the 
Tokai-to-Kamioka-and-Korea (T2KK) experiment, can observe the oscillation maximum 
of the i>n Ue transition probability at two vastly different base-line lengths. Our most 
important observation for the T2KK proposal is that the matter effect term in eq. (]14bp 
at Korea is about 3 times as large as that at Kamioka. Because the sign of the matter 
effect in A'^ follows the sign of A13, the amplitude of the z/^ Vf. oscillation for the normal 
hierarchy at IA13I ~ vr is larger than that for the inverted hierarchy. Because the matter 
effect is stronger at Korea, if the probability is found larger (smaller) at Korea than the 
one at Kamioka, we can conclude that the hierarchy is normal (inverted), irrespective 
of the other model parameters such as sin^ 26'rct and 5mns [14] . Let us examine this 
observation semi- quantitatively by using the approximate formulae in eqs. f|T3|) - f|T^ . 
The difference between the z/^ — > z/g oscillation amplitude at a far detector and that at 
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Figure 1: (a). The flux of tlie T2K z/^ beam for 10^^ POT/yr at L = 295 km for various 
off-axis angles between 0° and 3°. (b). The cross section per nucleon of the Charged 
Current Quasi Elastic (CCQE) events for z/g and z/^ off the water target, (c). The flux at 
295 km times z/g CCQE cross section, (d). Probability of z/^ Vf. transition at 295 km 
(solid line) and that at 1000 km (dashed line) calculated for m\ — m\ = 2.5 x 10~'^eV , 
ml - mj = 8.2 x 10-^eV^ sin^ 20atm = 1-0, sin^ 2^sol = 0.83, sin^ 2^rct = 0.10, 
i^MNS = 0°, and p = 2.8g/cm^ for L = 295 km, and p = 3.0g/cm^ for L = 1000 km. 
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(a) 




Figure 2: The off-axis angle of the neutrino beam from J- PARC on the sea level, when 
the beam center is 2.5° (a) and 3.0° (b) off at the SK site. The baseline lengths are 
shown by vertical contours, and the off-axis angles at the sea level are shown by elliptic 
contours between 0.5° and 3.0°. The SK is slightly off the corresponding contour because 
it is about 320 m above sea level. 
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a near detector is 



The error of the AP^ ^i^^(Ai3) can be estimated as 



(16) 



S ( AP,_,^ (A 



13, 



SPu^^uAL. 



+ 



5Pu^,^Ue{Lia.r) 




(17) 



Here is the number of z/g appearance events. The ratio between A^^'^'' and A''"'^'^'' at 
the maximum value of the oscillation probability, A13 = +7t, can be expressed as 



N- 



far 



Vfar ^'far (E^ at A13 = TT, Lfer) 



^^ear "^near 

{E^ at Ai3 = vr, L 

r 



(18) 



where V denotes the fiducial volume of the detector and ^{E^,, L) is the neutrino beam 
flux at L, which is proportional to (1/L)^. The cross section ratio at different energies 
drops out, because the neutrino cross section of CCQE events is almost constant in the 
0.7-5 GeV region; see Fig. [I](b). We therefore need to estimate the number of Ue CCQE 
events near the oscillation maximum at SK, A^near- We show in Fig. [3], typical numbers 
of expected CCQE events for the /i events (a) and the e events (b), for the 3.0° OAB 
at SK. The open squares show the expected numbers of events in a 200 MeV wide E^ 
bin, after 5 years (5 x 10^^ POT), at sin^ 26'rct = 0.1 and 5mns = 0° for the normal 
hierarchy, just as in Fig. [1] (d). From the two bins around Ei, = 0.6 GeV in Fig. [3]^b), 
we may estimate A^near ~ 130. 

The difference between the maximum value of the oscillation probability, APy^^^^, 
for the normal hierarchy (A13 = +7r) and that for the inverted hierarchy (A13 = — tt) 
can be expressed as 



AR 



AP. 



0.01 (t^^2^\ (j^ - ; 

V 0.10 ; V295km 



where we set Lnear = 295 km and used the approximation of eqs. ( |T3i) and f ll4b[) . The 
difference grows linearly with the distance, Lfar, as long as the oscillation maximum 
is covered by the flux. The significance of excluding the fake hierarchy can then be 
estimated as 



AP,^_,^ (+7r) - AP,^^,^ (-tt) 
5(AP,^_,,(+7r)) 

^sin^26'RCTV^V ^far 



2.3 



0.10 ; V295km 



- 1 



1 + 0.225 



lOOkt 



V 295km/ Vf 



far 



-1/2 



• (19) 
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Figure 3: The typical numbers of the /i events (a), and those of the e events (b), for the 
exposure time of 5 years (5 x 10^^ POT), for the 3.0° OAB at SK (open square), and 
for the 0.5° OAB at L = 1000 km with a 100 kt water Cerenkov detector (sohd circles). 
The input parameters are the same as those of Fig. [1] (d) [14]. 



We find that when we put a 100 kt detector at L = 1000 km, the significance can exceed 
the 3-0" level in this very rough estimate, which is confirmed in the following numerical 
studies. 

The phase-shift factor B*^ of eq. fllSp also contributes to the determination of the 
neutrino mass hierarchy pattern. Since sin^(x) is an even function of x, the magni- 
tude of IA13/2 -|- B'^l controls the —>■ oscillation in eq. f|T5l) . Because the value of 
IA13/2 + 5^1 varies by changing the sign of A13, we have the possibility to distinguish 
the sign of A13 by measuring the — > z/g oscillation phase. At Kamioka, however, 
is governed by the term proportional to cos(5mns; see eq. (I14d|) . and hence the fake hier- 
archy can reproduce the same IA13/2 -|- by changing the sign of cos^mns- In Korea, 
the magnitude of the matter effect term in eq. (114d[) is larger than that of the cos ^mns 
terms, and hence the fake hierarchy cannot reproduce the same oscillation phase. The ef- 
ficiency of these phase-shift contribution grows with the magnitude of 5^, which grows as 
cos Smns decreases. We find that the effect of the phase difference becomes as important 
as that of the amplitude difference at cos^mns ~ — 1- 

From the above consideration, we observe that it is useful to have a far detector in 
Korea in the region where the neutrino flux is significant around the z/^ — > Ug oscillation 
maximum, which is typically between 1.4 GeV to 3 GeV, see Fig. [IJd). Figs. [1] (a) and 
(c) show that beams at an off-axis angle smaller than 1° have this property. Figs. [21(a) 
and (b) show that such an off-axis beam will appear in a specific region near the east 
coast of Korea during the T2K experimental period. 

Once the neutrino mass hierarchy is determined, the T2KK experiment also has 
the ability to measure the leptonic CP phase uniquely without using the anti-neutrino 
beam. From eqs. (fT3l) and (114bl) . there are 2 unmeasured parameters which control 
the amplitude of oscillation, sin^6'RCT and sin^MNS- Due to the significantly 

different matter effect at IA13I = vr between Kamioka and Korea, we can constrain both 
sin^ 26'rct and 5mns uniquely from the two amplitudes. The value of cos 5mns is measured 
through the energy dependence of the z/^ v^. oscillation probability, through the phase 
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shift -B^; eq. (114dp . Since both sin^MNS and cos(5mns can be measured independently, 
the CP phase 5mns can be constrained uniquely. 



4 Analysis method 

Before we present the results of our numerical calculation, we would like to explain 
our treatment of the signals and background. In our case study, we consider a 100 kt 
level detector, in order to compensate for the decrease in the neutrino flux which is about 
(300 km/1, 000 km)^ ~ 1/10 of that at SK. We adopt a Water Cerenkov detector because 
it allows us to distinguish clearly the events from /x^ events. We use the CCQE events 
in our analysis, because it allows us to kinetically reconstruct the neutrino energy event 
by event. Since the Fermi- motion of the target nucleon dominates the uncertainty of 
the neutrino energy reconstruction, which is about 80 MeV, we take the width of the 
energy bin as SE^, = 200 MeV for E,y > 400 MeV, in the following analysis 0. The event 
numbers of from Ua flux {^u^) which is delivered by J- PARC [33] in the i-th energy 
bin, Ei = 200MeV x {i + 1) < E^ < E^ + 5Ey, are then calculated as 

= MNa j^^ {E) P,„^,^ {E) a«^(E) dE , (20) 

where Py^^y^ is the neutrino oscillation probability including the matter effect, M is the 
detector mass, Na = 6.017 x 10^^ is the Avogadro constant, and a^^ is the CCQE cross 
section per nucleon in water [4]. All the primary as well as secondary fluxes used in our 
analysis are obtained from the website [33]. The fiducial volume of Super-Kamiokande 
is 22.5 kt, and we assume that a detector in Korea is 100 kt. CCQE events have been 
selected as single-ring events at K2K. Efficiency of the technique has been estimated 
to be about 80% around 600 MeV and about 60% around a few GeV [10,34]. In the 
following numerical analysis we set the efficiency to be 100% at both SK and at a far 
detector for brevity [§. 

When the secondary neutrino flux (z/g , Ve , ^p) of the z/^ primary beam is considered, 
the e-like and //-like events for the i-th bin are obtained as 

K = Kii^,)+ E ^:M, (« = e,/i). (21) 

The second term in eq. (!2T]) corresponds to the contribution from the secondary neutrino 

^At a few GeV region, contributions from soft-charged tt emission processes to the CCQE signal 
events become significant [34]. Those events have soft charged tt's which do not emit Cerenkov hghts, 
and hence the reconstructed neutrino energy is underestimated by about 300 MeV [10]. In this analysis 
we do not consider contribution from non-CCQE events, and leave the studies of their impacts for the 
future. 

^ Our results can be regarded as those for larger fiducial volumes, such as 28.1 kt at SK and 167kt 
at a far detector, if we take the efficiencies of 80% and 60% estimated in refs [10,34]. 
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flux. In our analysis, we calculate NJ^^^ by assuming the following input parameters: 



\ml - m2p°P"t = 2.5 x lO'^ eV^ 
(ml - m'^Y'^P''' = 8.2 x 10"^ eV^ 



Sm P^TM 



sm' 



26, 



input 
SOL 



0.5, 
= 0.83, 



(22) 



with the constant matter density, p™?*^* = 2.8 g/cm^ along T2K and p™P"* = 3.0 g/cm^ 
for the Tokai-to-Korea baseline, which goes through deeper in the earth than that of 
T2kI2|. We use the above numerical values in order to keep the consistency with our 
previous report [14]. Our results are not sensitive to the small change in the input values 
in eq. fl22|) but the quantitative results on the mass hierarchy determination depend on 
the mater density along the baseline. Dedicated study on the impacts of the matter 
density profile along the T2KK baseline will be reported elsewhere. We examine various 
input values of sin^26'RCT , ^mns ? and the sign of m| — Since our concern is the 
possibility to distinguish the neutrino mass hierarchy and to constrain the CP phase 
uniquely, we study how the above 'data' can constrain the model parameters by using 
the function 

^X^ = xIk + Xir + Xsys + Xpara • (23) 

Here the first two terms, Xsk ^^nd Xkt^ measure the parameter dependence of the fit to 
the SK and the Korean detector data. 



XsK.Kr — 



(24) 



a=e,fi 



where the summation is over all bins from 0.4 GeV to 5.0 GeV for A^^, 0.4 GeV to 1.2 
GeV for A^'^ at SK, and 0.4 GeV to 2.8GeV for at Korea. Here N^^ is the calculated 
number of events in the i-th bin, and its square root gives the statistical error. We 
include the contribution of the /i-events in order to constrain the absolute value of A13 
strongly in this analysis, because a small error of A13 dilutes the phase shift [14,25]. 

A^f* is calculated by allowing the model parameters to vary freely and by allowing 
for systematic errors. In our analysis, we consider 4 types of systematic errors. The 
first ones are for the overall normalization of each neutrino flux, for which we assign 3% 
errors, 

= 1±0.03, (25) 

for all neutrino flavor, which are taken common for T2K and the Tokai-to-Korea experi- 
ment. The second systematic error is for the uncertainty in the matter density, for which 
we allow 3% overall uncertainty along the baseline, independently for T2K (/p^) and 
the Tokai-to-Korea experiment {f^"^) ■ 



pf = flpr' (^ = SK, Kr). (26) 

•^Because the sea between Japan and Korea is less than 1km deep, there is no contribution from the 
water to the average density. 
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The third uncertainty is for the CCQE cross section. Since z/g and z/^ CCQE cross 
sections are expected to be very similar theoretically, we assign a common overall error 
of 3% for z/g and z/^ {fe^ = fj}^ = ff^), and an independent 3% error for z/g and 



CCQE cross sections {fi^ = f^^ = f^)- The last one is the uncertainty of the fiducial 
volume, for which we assign 3% error independently for T2K (/v*^) and the Tokai-to- 
Korea experiment (/vO- N^^^ is then calculated as 



and Xsys has four terms; 

2 / nn-p \ 2 




""'^^ .=£,.(^0-03 ) ^|r,r(^"0.03 ) 

(28) 

In short, we assign 3% errors for the normalization of each neutrino flux, the z/; and 
ui CCQE cross sections, the effective matter density along each base line, and for the 
fiducial volume of SK, and that of the Korean detector. 

Finally, Xpara accounts for the external constraints on the model parameters: 



^P"^'^ 0.6 X 10-5 J ^[ 0.07 J 

^ |^ sin^2gfcT-sin^2<g^y _ ^^g) 

The first two terms correspond to the present experimental constraints summarized in 
eq. In the second term, we adapt the larger error than that shown in eq. (15bp 

for continuity of the previous work, ref. [14]. In the last term, we assume that the 
planned future reactor experiments [26,27] should measure sin^ 26'rct with the expected 
uncertainty of 0.01, during the T2KK experimental period. In total, our Ax^ function 
depends on 16 parameters, the 6 model parameters and the 10 normalization factors. 

5 Determination of the mass hierarchy pattern 

In this section we show the results of our numerical calculation. First, we search for 
the best combination of the off-axis angle at SK and that at a Korean detector, as well 
as the base-line length up to Korean detector for determining the sign of ml — ml. 
For this purpose, we first calculate the expected number of the z/^ z/g CCQE events 
at both detectors by assuming either normal or inverted hierarchy, and then examine 
if the resulting 'data' can be fitted for the opposite hierarchy by adjusting the model 
parameters. 
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Figure 4: Minimum Ax^ of the T2KK two detector experiment after 5 years of running 
(5 X 10^^ POT) as functions of the off-axis angle and the base-hne length of the far- 
detector from J-PARC at Tokai, when the normal hierarchy (m| — mf > 0) is assumed 
in generating the events, and the inverted hierarchy (mg — mf < 0) is assumed in the 
fit. The left-hand figure (a) is for the 2.5° OAB at SK, and the right-hand one (b) is for 
the 3.0° OAB beam at SK. The input parameters are the same as those of Fig. [T] (d); in 
particular, sin^ = 0-10 (J^ns — 0°- 

We show in Fig. H] the minimum Ax^ expected at the T2KK two detector experiment 
after 5 years of running (5 x lO^^POT), as functions of the off-axis angle and the base- 
line length of the far-detector site from J-PARC at Tokai, when the normal hierarchy 
{ml— ml > 0) is assumed in generating the events, and the inverted hierarchy {m\—m\ < 
0) is assumed in the fit. The left-hand figure (a) is for the 2.5° OAB at SK, and the 
right-hand one (b) is for the 3.0° OAB beam at SK. The input parameters are choose as 
in eq. (12^ . sm^29^^Q^ = 0.10 and ^J^^s ~ '^^^ ^^^^ symbols, solid circle, open circle, 
triangle, and square are for L = 1000km, 1050km, 1100km, and 1150km, respectively. 
There are no data points at 0.5° in Fig. 111(a) for the 2.5° OAB at SK, because the 0.5° 
off-axis beam does not reach Korea: see Fig. [2]^a). It is clearly seen from Fig. H] that 
the best combination of off-axis angles are 3° for SK and 0.5° for the Korean detector at 
L = 1000 km. The 0.5° off-axis beam has strong flux up to ~ 2.2 GeV, which overlaps 
significantly with the broad peak of the z/^ v^. oscillation at L = 1000 km; see Fig. [1] 
(a), (c) and (d). Because the number of the Vf. CCQE events is large enough around the 
oscillation maximum for sin^ 2^rct ~ 0.1, both at SK and at the far detector in Korea, 
we are able to measure the difference in the magnitude of the z/^ — z/e probability at 
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Figure 5: The same as Fig. HI but when the input data is calculated for the inverted 
hierarchy (m| — ml < 0) and the fit is performed by assuming the normal hierarchy 
(m| — ml > 0). All the other parameters are the same as those in Fig. HI 

two vastly different baselines, and can hence distinguish between the normal hierarchy 
and the inverted hierarchy. We can reject the fake hierarchy at 4.7-0" level in our simple 
simulation with this combination of 3.0° at SK and 0.5° at L = 1000 km. If we remove the 
the constraint of the future reactor experiment, the last term in eq. (I29p . the minimum 
Ax^ value drops from 23 to 18, for the combination of 3.0° at SK and 0.5° at L = 1000 
km; see [14]. 

We find from Fig. |4] that the 1.0° OAB in Korea still keeps the sensitivity to the 
neutrino mass hierarchy, where both a combination of 1.0° at L = 1000 km and 2.5° 
OAB at SK (Fig. 1 (a)) and that of 1.0° at L = 1000 km ~ 1150 km and the 3.0° OAB 
at SK (Fig. m (b)) distinguish the neutrino mass hierarchy nearly at 4-0" level in our 
simulation. This is because the CCQE cross section times the flux of 1.0° OAB extends 
to ~ 1.7 GeV, see Fig. [1] (a) and (c), which barely overlaps with the broad peak region 
of the z/^ — i> z/e oscillation probability shown in Fig. [I](d). From Fig. [2] (a) and (b), we 
find that the 1.0° OAB is observable only in the east coast of Korea {L ~ 1000 km) for 
the 2.5° OAB at SK, whereas for the 3.0° OAB at SK, it can be observed at various 
base-line lengths up to ~ 1150 km. The small values of Ax^ for larger off- axis angles in 
Fig. m tell us that it is essential to choose the location of the detector in Korea where 
the off-axis angle is smaller than 1.0°. 

Fig. O is the same as Fig. HJ but when the input data is calculated for the inverted 
hierarchy and the fit is performed by assuming the normal hierarchy. It is remarkable 
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Figure 6: Capability of the T2KK two-detector experiment to determine the neutrino 
mass hierarchy, (a) when the mass hierarchy is normal {m\ — m\ > 0), and (b) when it is 
inverted {ml—ml < 0). The numerical results are obtained for a combination of 3.0°OAB 
at SK and 0.5°off-axis at L = 1000km with a 100 kt water Cerenkov detector, after 5 
years of running (5 x 10^^ POT). In each figure the input data are calculated for the model 
parameters at various sin^ 26*^^* and ^J^ns ; been performed by surveying 

the whole parameter space with the opposite mass hierarchy. The resulting values of 
minimum Ax^ are shown as contours for 2, 3, 4, and 5a. The wrong hierarchy can be 
excluded with the corresponding confidence level if the true sin^ 26'rct and (5mns values lie 
in the right-hand side of each contour. The model parameters are set at (m| — mf)'°P"* = 
2.5 X 10-^2 (a), -2.5x lO-^eV^ (b), (m^ - m?)*''P'^* = 8.2 x 10-5eV^ sm^2e'^^^ = 1.0, 
sm^29fSl^ = 0.83, p'°P"* = 2.8g/cm^ for SK, and p'^^P"* = 3.0g/cm^ for L = 1000km. 

that almost the same level of the capability to distinguish the neutrino mass hierarchy 
can be achieved for the OAB ^ 1° at L ~ 1000 km even when the hierarchy is inverted 
— ml < 0). Slight decreases of the minimum Ax^ value of the best combinations, 
23.5 in Fig. Mh) to 21.3 in Fig. EJ^b), and 14.8 in Fig. Ka) to 14.6 in Fig. EJ^a), can 
be attributed to the smaller expected number of the z/g appearance events because of 
the matter effect which suppresses the probability in the inverted hierarchy. We may 
conclude that a far detector that observes the T2K neutrino beam at L ~ 1000 km and 
the off-axis angle ^1.0° has the potential to determine the neutrino mass hierarchy 
whether it is normal or inverted. 

Because we find from Fig. HJand Fig. [5] that a combination of 0.5° OAB at L ~ 1000 
km and the 3.0° OAB at SK has a significantly better capability of determining the neu- 
trino mass hierarchy, we study in the following physics potential of this preferred T2KK 
set up in more detail. In particular, we investigate the whole un-explored parameter 
space of the three neutrino model. 

First in Fig. El we summarize our findings on the capability of the T2KK exper- 
iment to determine the neutrino mass hierarchy in the whole space of sin^ 2^rct and 
5mns • Fig. [6](a) shows our result when the mass hierarchy is normal (mg — ml > 0), 
and (b) when it is inverted (mg — mf < 0). In each figure the input data are calcu- 
lated for the model parameters at various sin^ 29^rct and 5mns points, and the fit has 
been performed by surveying the whole parameter space, but under the opposite mass 
hierarchy. The resulting values of the minimum Ax^ are shown as contours for 2, 3, 4, 
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and 5-cr. The wrong mass hierarchy can be excluded with the corresponding confidence 
level if the true sin^ 26'rct value lies in the right-hand side of each contour along the 
true value of 5mns ("^mns)- particular, the minimum Ax^ values of 22 for the point 
(sin^ '^MNs) ~ (0.10, 0°) in Fig. [H](a) corresponds to the highest point in Fig.Hl^b), 

and the corresponding value of 21 in Fig. [6](b) is the highest point in Fig. [5]^b). We find 
that the wrong hierarchy can be excluded at the 3-a level if sin 

2^RCT > 0.05 (0.06) if 

the hierarchy is normal (inverted). 

It is remarkable that the ^mns = 0° case chosen to plot Fig. [Hand Fig. [S] turns out to 
be the case when it is most difficult to determine the neutrino mass hierarchy. If Suns = 
180°, the wrong hierarchy can be excluded at the 3-a level for sin 

2<CT ~ 0.02 for the 

normal hierarchy (Fig. [6t^a)) or sin^ 26'^^^* ~ 0.03 for the inverted hierarchy (Fig. El^b)). 
The origin of the 5mns dependence is the difference of the oscillation phase at the far 
detector in Korea. From eqs. fll4dp and f[T^ . the difference of the oscillation phase near 
the oscillation maximum, |Ai3| = n, between the input and the fit is expressed as. 



^■^ _|_ Re 
2 ' input 



A 



13 



B 



fit 



±0.15 < cos(5J^^g* 



0.10 y/^ 

- 1 + cos 



fit 

MNS 




^0.58 ( — ^ . 

VlOOOkmy 



(30) 



The upper sign is for the normal hierarchy, and the lower sign is for the inverted hierarchy. 
The phase-shift difference depends on both cos(5J^ns cos(5^iv^g. As explained in 
section 3, below eq. (IT^ . when cos^J^^s* ~ ^ ('^mns ~ 0°) the phase shift is smaller 
than that with the other ^J^ns at L = 1000 km. Therefore the fitted value of cos^mns 
(cos5^Ns) ^Is*^ tends to be large and has the opposite sign for the fake hierarchy. If 
cos5j^NS* ~ —1, it is not possible to compensate for the phase-shift difference of eq. fl30l) 
even by choosing cos^^^s ~ ^^"^ the significantly higher minimum Ax^ value results 
in Fig. [H](a) and (b). In general, cos^^j^fg > is favored even when cos^JJJ^s* < in order 
to minimize the phase-shift difference of eq. ( 1301) . This is clearly seen in Fig. [71 where we 
show the values of the best fit parameters, sin^ '^(^rct ^mns; ^t the minimum Ax^ 
point of the analysis of Fig. [61 The results for the normal hierarchy (a) and those for the 
inverted hierarchy (b), are shown correspondingly to the fit of Fig. M^a) and Fig. Mjo), 
respectively. The thick vertical lines are the sin 

2e|*CT contours at 0.02, 0.06, 0.10 and 
0.14. The thinner contours give the S^-^^ values. We find that the value of 5^ns around 
0° is almost always favored as expected. 

Here let us try to explain more detailed features of Fig. [6] and Fig. [7] by separating 
the parameter space of sin^ 26*1^^* and 5mns i^to 4 regions. 

1. small sin^ 2^^CT region ( sin^ 26*1^^^* < 0.04) at any ^MNg: In this region the 
phase difference eq. f[30|) is mainly controlled by the cos 5mns terms, because of the 

1 / J sin^ 26'rct enhancement over the matter effect term. It is hence relatively easy 
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(a) input : normal, fit : inverted (b) input : inverted, fit : normal 




Figure 7: The values of the fit parameters, sin^26'^Qrp and 5mns' minimum Ax^ 

point of the analysis of Fig. [6] are shown. The results for the normal hierarchy (a) and 
those for the inverted hierarchy (b) are shown correspondingly to the fit of Fig. [6] (a) 
and Fig. [6] (b), respectively. The thick vertical lines are the sin^ 

2^fcT contours at 0.02, 

0.06, 0.10, and 0.14. The thinner contours give the 5mns values. 

to make the difference small by adjusting cos ^^^g + cos (5J^ns ~ 0- The hierarchy is 
determined essentially by the difference of the z/^ v^. oscillation amplitude only. 

2. sin^26'j^CT ~ 0-04 at ^mns ~ 180°: Although the effect of cos^J^^s ~ ~^ can- 
celed by choosing cos^^^s ~ +-'-5 difference from the matter effect term in 
eq. (!30l) cannot be canceled. Therefore in this region the hierarchy is determined 
by the differences of both the amplitude and the oscillation phase. 

6. sm 

2^RCT ~ 0.04 at 

"^MNS ~ ±90°: In eq. (ED]), the difference is controlled by the 
matter effect term and the cos5^j,^g term because cos^J^^s ~ 0- ^^i^ region, we 
can make the phase-shift difference small by choosing cos^^^g > 0. 

4. sin^26'J^CT ~ 0-04 to 0.08 at 5mns ~ 0°- ^^is region the phase-shift difference 
eq. ( l30l) at IA13I ~ vr can be made small at cos(5mns ~ ^1 but the difference at 
IA13I ~ 27r becomes large. Because the flux of 0.5° off-axis beam is strong at lower 
energies where tt < IA13I < 27r, the growth of the phase-shift difference eq. (!30l) at 
larger IA13I cannot be compensated. This explains why the minimum Ax^ value 
in this region is larger than the one for the case 3. 

The systematics of the oscillation phase is rather complicated, but its effect turns out to 
be significant in determining the neutrino mass hierarchy. 

Before closing the section, let us briefly study the value of sin^ 26'{^Qrp in Fig. [71 In 
Fig. [Zl(a), sin^ 2^|*Qrp is larger than sin^ 26'^^^ 5 whereas in Fig. [7](b), sin^ 2^|*Qrp is smaller 
than sin^2^|^CT- This is because the same oscillation amplitude can be obtained by 
choosing sin^ '^^bIgt: > sin^ '^dRcr when the hierarchy is normal but it is assumed to be 
inverted in the fit, and vice versa for the opposite case. The fitted value of sin^ 26*^^ 
cannot deviate too much from the input value of sin^ 26*^^5 however, because of the 
constraint from the proposed reactor experiments, according to the last term in eq. fl29|) . 
and also because of the SK measurement of the z/^ —>■ Ve probability, which is much less 
sensitive to the mass hierarchy difference. 
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6 Measurement of the CP phase 



In this section we investigate the measurement of sin^ 2^rct and 5mns for our preferred 
combination of the 3.0° OAB at SK and 0.5° OAB at L = 1000 km. This combination of 
the T2KK experiment allows us to measure the oscillation around the oscillation 

maximum at two base-line lengths, which can be parametrized as in eq. ( JT3l) . in terms 
of the amplitude shift eq. (114bp and the phase shift eq. (114dp . Once the neutrino mass 
hierarchy is determined as explained in the previous section, the terms proportional to 
|Ai3| /tt in the amplitude shift eq. fll4bl) measure sin^MNS, and those in the phase shift 
eq. fll4dl) measure cos^mns- In the T2KK two detector system, both sin^ 26'rct and 
sin^MNS can be determined uniquely because the amplitude shift eq. fll4b|) has signifi- 



cantly different matter effect contributions between SK and the far detector. The phase 
shift measurement of the term eq. fll4dl) constrains cos(5mns independent of sin^MNS- 

We show in Fig. [8] and Fig. [9] regions allowed by this experiment in the plane of 
sin^ 26'rct and (5mns- The mean values of the input data are calculated for the parameters 
of eq. (j22l) . In each figure, input points (sin^ 2^^^^ 5 '^mns ) shown by sohd-circles for 
sin^ 26'J^qJ^* between 0.02 and 0.1, with an interval of 0.02, and for four values of ^mnsj 
0°(a), 90° (b), 180° (c), and — 90°(d). The regions where the minimum Ax^ value is less 
than 1, 4, 9 are depicted by solid, dashed, and dotted boundaries, respectively. Fig. [8] 
is for the normal hierarchy, and Fig. [9] is for the inverted hierarchy. From these figures, 
we find that Sum can be constrained to ±30° at 1-cr level, when sin^ 26*^^ ~ 0-02 as 
long as the neutrino mass hierarchy is determined. If we remove the last term in eq. fl2^ . 
the error of 5mns changes little but that of sin^ 26'rct grows significantly; see Fig. 5 of 
Ref. [14] 

As shown in Fig. [6](a) and (b), the mass hierarchy cannot be determined at 3-0" level 
(Ax^ > 9) when sin^ 2^^^^ is too small. In case of the input parameters of Fig. [8] 
for the normal hierarchy, this is the case for sin^ 26'^^ = 0-02 at (5mns ~ 0° l^); 
sin^ 26'^^^ = 0.02 and 0.04 at 5^ = 90° (b), and sin^ 2^1^;^:^,* = 0.02 at 5'^^^ = -90°(d). 
For those input points, there appear an additional allowed region whose center (local 
minimum of Ax^) is shown by a solid square. No extra allowed region appears for 
(^MNS = 180° in Fig. [8](c), in accordance with the result of Fig. [6](a). In case of Fig. [9] for 
the inverted hierarchy, the local minimum appears for sill ^^p^^rp — 0.02 cit (^jvi^s — ^ 
(a), siYi^2e^l^ = 0.02 and 0.04 at 5^ = 90° (b), sin^ 2^|^p^* = 0.02 at 5^ = 180° 
(c), and sin^ 2^|^p^* = 0.02 at 5^ = -90°(d). 

It is remarkable that the error of (5mns is almost independent of sin^ 20^^* value 
between 0.02 and 0.1, for all the four input values of ^mnS) 0°, ±90°, and 180°. This 
is remarkable because the event number iVg is proportional to sin^ 2^rct according to 
eq. fll3p . and hence the statistical error of the measurement of the amplitude and the 



phase should be proportional to or l/ysin 6'rct- This increase in the error for 

small sin^ 2^^^^* values is canceled by the increased sensitivities of both the amplitude 
and the phase shift to sin^MNS and cos^mns? respectively, which are both proportional 



to l/ysin^ ^RCT- The two effects cancel rather accurately, and we find that the error of 
(^MNS is almost independent of the input values of sin^ 2^rct and (5mns- 
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Figure 8: Capability of the T2KK two detector experiment for measuring sin^ 2^rct and 
(^MNS ■ Allowed regions in the plane of sin^ 2^rct and 5mns are shown for a combination 
of 3.0°OAB at SK and 0.5°at L = 1000km with a 100 kt water Cerenkov detector after 
5 years of running (5 x lO^^ POT). The input values of sin^ 2^rct are 0.02, 0.04, 0.06, 
0.08 and 0.10 for (5mns = 0°(a), 90°(b), 180°(c), and -90°(d). The normal hierarchy is 
assumed at m| — mf = 2.5 x 10~^eV^, and the other parameters are the same as those in 
Fig. O The input points are shown as solid blobs, where Ax^ = by definition. The 1-, 
2-, and 3-0" contours are then shown by solid, dashed, and dotted hues, respectively. For 
the input values of {sin^ 2e'^l^\ ) = (0-02, 0°) (a), (0.02, 90°) and (0.04, 90°) (b) 
and (0.02, —90°) (d), there appear additional allowed regions when the mass hierarchy 
is chosen with the wrong sign in the fit, where the local minimal Ax^ point is depicted 
by a solid square. 



22 




90 



-90 



00 0.02 


0.04 


0.06 0.08 

sin^20RCT 


0.10 


0.12 


0.14 




■'■ ~> 






















(6) ■ 






















































































































































































































































































3 










^5 













0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 

sin220RCT 



-180 I ' ' ' ' ' ' ' ' ' ' ' ' ' ^ 

0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 

sin^2eRCT 





























(CI) 




















































































































































































































1 


... .. . 


i 
















i 


m 





































































0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 

sin220RCT 



Figure 9: The same as Fig. [HI but when the events are calculated for the inverted 
hierarchy, i.e., m| — mf = —2.5 x 10~'^eV^. Just like in Fig. [HI additional allowed 
regions, when the wrong sign of the ml — ml is chosen in the fit, appear for all the ^mns 
cases at sin^ 29'^^^ = 0.02, and for = 90° (b) at sin^ 2^1;^^^ = 0-04. 

7 Conclusion and discussions 

In this paper we study physics potential of the T2KK proposal [12-14], where a far 
detector along the T2K neutrino beam line is placed in Korea. We find that the off-axis 
neutrino beam from J-PARC at Tokai village for the T2K project has significant intensity 
at a few GeV range when the far detector is placed in the east coast of Korea where 
the beam at less than 1.0° off-axis can be observed at a base-line length of L ~ 1000 
km. The resulting two detector system can observe the Ve oscillation probability 

near the oscillation maximum at two different energies, if sin^ 2^rct is not too small. 
We examine, in particular, the capability of determining the mass hierarchy pattern and 
the CP phase of the lepton-flavor- mixing matrix when a 100 kt water Cerenkov detector 
is placed at various locations in Korea for the off-axis beam (OAB) of 2.5° and 3.0° 
at the Super-Kamiokande site. The best results are found for a combination of 3.0° 
OAB at SK (L = 295km) and 0.5° OAB at L = 1000km, where the mass hierarchy 
pattern can be determined at 3-0" level for sin^ 20rct ~ 0.05 (0.06) when the hierarchy is 
normal (inverted), after 5 years of running (5 x 10^^ POT). The sensitivity of the T2KK 
experiment on the neutrino mass hierarchy depends not only on sin^ 2^rct but also on 
^MNS • We explore the sensitivity in the whole space of sin^ 2^r,ct and ^mns , and the 
results are shown in Fig. [6]^a) for the normal hierarchy and in Fig. [6]^b) for the inverted 
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hierarchy. Significantly higher sensitivity is found for Suns ~ 180° for both hierarchy 
cases. 

We also find that the leptonic CP phase, 5mns 5 can be constrained uniquely, without 
invoking anti-neutrino beams, as long as the mass hierarchy pattern is determined: see 
Fig. [8] for the normal hierarchy, and Fig. [9] for the inverted hierarchy. 

Those results are obtained by assuming that the neutrino energy can be reconstructed 
with a hundred MeV uncertainty for the charged current quasi-elastic events, and the 
earth matter density along the baseline can be determined with 3% accuracy. All our 
numerical results have been understood semi-quantitatively, by using the approximative 
expression for the and z/^ — > z/g oscillation probabilities, eqs. ffTOj) . ffTTl) . f|T3l) 

f|T^ . and f|T5|) . where only the linear and quadratic terms in the matter effect and those 
in sin^ 26'rct are kept in the oscillation probabilities. 

Our results are based upon a very simple treatment of the systematic errors where 3% 
overall errors are assigned for all the 10 normalization factors of eq. fl25]) . We find that the 
significance of the mass hierarchy determination is not affected much even if we enlarge 
all the systematic errors to 10% except for the matter density uncertainties. This means 
that the errors of the T2KK experiment proposal in this paper are dominated by the 
statistical error. Therefore, if we make the detectors at Kamioka and/or Korea larger, or 
if the intensity of the J-PARC beam is stronger, the significance of the measurement will 
grow as proportional to the square root of the product of the volume, the intensity, and 
the exposure time. It is not so clear, however, what is the best volume ratio between 
the detector in Kamioka and that in Korea. We show in Fig. [10] the minimum Ax^ 
as functions of the volume ratio of the near (Kamioka) and far (Korea) detectors while 
keeping the total volume at 600 kt for 5 x 10^^ POT. We assume the normal hierarchy for 
the input and the inverted hierarchy in the fit. We examine 8 cases, for sin^ 2^^^^^* = 0.1 
(solid lines) and sm^2e''^^^ = 0.06 (dotted hues), and for 5^ = 0° (a), 90° (b), 180° 
(c), and — 90°(d). It is clearly seen from Fig. [10] that a 600 kt detector in Kamioka alone 
cannot resolve the mass hierarchy at all, because there is a little difference in z/^ — >■ z/g 
transition probability between the normal hierarchy and the inverted hierarchy. On the 
other hand, in the case of only a Korean detector with 600 kt, the minimum Ax^ value is 
not much smaller than the best case. This is because the constraint of sin^ 2^rct from the 
future reactor neutrino experiment replaces the role of the near detector which measures 
the Ufj^ z/g transition at low energies where the matter effect is small. We find that the 
minimum Ax^ value of 23.5 in Fig. [TI]r a) at the volume ratio of 1:5 (~ 22.5 : 100) is about 
4.5 times as large as the minimum Ax^ value in Fig. Il](b), confirming the dominance 
of the statistical error in our analysis. If we request that the minimum Ax^ should be 
at least 80% of its optimal value, then the near-to-far volume ratio should be between 
0.5 : 5.5 and 2.5 : 3.5. More volumes should be given to the far detector than to the 
near detector. The above results are obtained with a 100% detection efficiency for the 
CCQE events. If we set 80% efficiency for SK and 60% efficiency for a far detector in 
Korea [34], the value of Ax^ in Fig. |l]for the combination of 3.0° for SK and 0.5° for 
a far detector placed at L = 1000 km drops from 22.6 to 14.5, roughly scaling with the 
effective number of events at a far detector. 

We note here that our proposal is significantly different from that of Ref. [13] in the 
following aspects. 

1. We propose to use a combination of large off- axis angle at SK and a small off-axis 
angle at the far detector, so that the first oscillation maximum can be observed at 
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Figure 10: Optimal ratio of the fiducial volumes of two detectors, one at Kamioka 
[L = 295 km) and the other at L = 1000 km in Korea, for determining the neutrino 
mass hierarchy. The input data are calculated for the normal hierarchy at mg — ml = 
2.5 X 10~^eV^ and the minimum Ax^ of the fit with the wrong hierarchy is shown for 
sm'^2e''^^^^ = 0.10 (sohd lines) and 0.06 (dotted lines) and for ^j^^f = 0° (a), 90° (b), 
180° (c), and 270° (d), when the sum of the fiducial volumes is fixed at 600 kt. The other 
parameters are same as those in Fig. [T](d), and the results are calculated for 5x lO^^POT. 

both detectors. 

2. As a consequence, the CP phase can be determined uniquely without invoking 
experiments with an anti-neutrino beam. 

3. We propose that the far detector should have larger fiducial volume in order to 
maximize the physics outputs, see Fig. [TUJ 

The use of two different off-axis angles implies that our proposal should suffer from 
larger systematic errors, especially for high energy events at the far detector. Among 
the potentially serious background which we could not estimate in this paper are; 

• possible miss-identification of NC 7r° production as CCQE events, 

• possible miss-identification of soft vr emission events as CCQE events. 

Although the above uncertainties were found to be rather small at K2K experiments [35], 
we should expect them to be more serious at high energies. Dedicated studies on the fake 
reconstructed neutrino energy of the NC tt" events and those on the correlation between 
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the true and the reconstructed neutrino energy of the CC soft tt events are mandatory. 
In addition, careful studies including possible energy dependence of the flux and cross 
section uncertainties, and also the location dependence of the matter density, may be 
needed to justify the physics case of the T2KK proposal. 
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